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Monte-Carlo  modeling  of  the  experimentally  determined  transition 
temperature^^ ^ as  a function  of  coverage  is  reported  for  the  chemisorp- 
tion system  W(110)p(2xl )-0.  This  system  is  one  of  the  first  to  be 
studied  in  detail  that  is  "closed";  i.e.,  it  is  in  equilibrium  neither 
with  the  gas  phase  nor  with  0 dissolved  into  bulk  W.  Thus  a determina- 
tion of  the  temperature-coverage  phase  diagram  presents  the  possibility 
of  studying  the  thermodynamics  of  two-dimensional  systems  and  through 
this  the  study  of  adatom-adatom  interactions  that  lead  to  the  formation 
of  two-dimensional  ordered  phases.  Additionally,  of  course,  through 
use  of  the  lever  rule  the  concentration  of  adsorbed  atoms  in  the  ordered 
vs  the  disordered  phase  at  any  temperature  and  coverage  can  easily  be 
determined. 

The  measurements  to  which  the  Monte  Carlo  calculations  are  com- 
pared were  taken  in  a simple  LEED  diffractometer  consisting  of  gonio- 
meter, Faraday  cup  collector,  and  electron  gun.  The  angular  distribution 
of  intensity  in  a superlattice  [e.g.  (1/2  1/2)]  reflection  was  measured 
as  a function  of  coverage  and  temperature.  The  peak  intensity  can  be 

related  to  the  transition  temperature,^^ ^ while  the  angular  width 

(21 

reflects  the  size  of  the  ordered  regions,  or  islands.'  ' 

The  ordering  of  overlayer  structures  in  general  involves  several 
interactions.  Thus  in  modeling  this  ordering,  a variety  of  experimental 
inputs  is  required.  In  the  past,  usually  only  the  transition  temperature 
at  fixed  coverage  was  fitted.  A better  approach  is  to  fit  transition 
temperatures  as  a function  of  coverage,  but  even  here,  a number  of 
models  with  different  interactions  are  satisfactory.  An  additional 
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experimental  quantity  that  can  be  fitted  is  the  island  shape  at  low 
coverage.  The  present  results  demonstrate  that  the  latter  is  a useful 
parameter  in  selecting  the  proper  interactions. 

Monte-Carlo  calculations  were  performed  for  several  models 
of  the  adatom-adatom  interactions  in  the  p(2xl)-0  overlayer  on  W(llO). 

The  first  corresponds  to  the  model  used  in  an  analytical  determination 
of  the  adatom-adatom  interactions. ^ The  others  represent  two  addi- 
tional possibilities  for  the  interactions.  All  gave  good  agreement 
with  the  transition  temperature  and  varying  but  still  reasonable 
agreement  with  the  decay  of  the  superlattice  beam  intensity  (proportional 
to  the  square  of  the  order  parameter).  All,  however,  gave  incorrect 
island  shapes. 

The  experimentally  observed  diffraction  features  are  round, 

implying  that  on  the  average  the  islands  are  round.  Since  only  two 

orientations  are  possible  by  symmetry,  and  these  do  not  interfere  with 

each  other,  this  implies  that  individual  islands  are  round.  Model  1 

(see  Fig.  1)  gave  islands  elongated  along  the  open  direction,  implying 

the  net  attractive  interaction  along  this  direction  was  chosen  relatively 

too  strong  in  the  model.  Model  2 gave  a slight  elongation  in  the 

close-packed  direction,  implying  relatively  too  much  net  attractive 

(3) 

interaction  along  this  direction.  Model  3 had  a poor  tendency  to  order 
at  temperatures  where  the  other  two  were  well-ordered,  and  indicated 
row  formation  as  in  Model  2. 

Figure  2 shows  the  calculated  island  size  distribution  at  two 
different  temperatures  and  a coverage  of  6 = 0.2  for  Model  2.  At  the 
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lower  temperature,  the  adsorbed  atoms  are  more  or  less  in  one  ordered 

island,  with  an  equilibrium  density  of  adatoms  in  the  "sea".  The 

tendency  toward  elongation  of  the  island  mentioned  above  is  evident. 

At  the  higher  temperature,  several  smaller  islands  exist,  with  both 

orientations  now  present.  The  size  of  these  islands  corresponds  quite 

(21 

well  with  the  experimentally  observed  sizes.'  ' 

The  results  as  listed  in  Figure  1 still  have  too  wide  a range 

of  uncertainty  for  meaningful  comparison  to  calculations  of  the  indirect 

(4  51 

oscillatory  interaction  in  chemisorbed  layers.'  ’ ' However,  in  this 

-5  (4  51 

particular  system  an  r decay'  ’ ' appears  to  be  too  rapid.  An  exten- 
sion of  the  experimentally  determined  phase  boundaries  to  lower  coverages, 
as  well  as  continued  Monte-Carlo  modeling  to  improve  the  fit  to  island 
shapes,  is  in  progress. 

Further  Monte-Carlo  calculations  are  also  underway  with  a goal  of 
providing  a qualitative  description  of  the  phase  boundaries  over  the 
full  coverage  range  0 1 0 1 . For  coverages  greater  than  0 = 0.5,  the  phase 

diagram  looks  quite  different^^^  than  for  0 < 0.5.  To  fit  these  data 
requires  the  inclusion  of  three-body  forces  in  the  calculations.  Detailed 
results  of  this  work  will  be  reported  elsewhere. 
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Figure  1 


Figure  2 


FIGURE  CAPTIONS 


Models  of  adatom  interactions  in  the  p(2xl)  0 overlayer. 
Model  1 is  the  same  as  the  analytical  model  of  Ref.  1, 
and  is  the  one  with  the  least  number  of  parameters  for 
a bridge-bonded  site.  Models  2 and  3 are  appropriate 
for  the  3 -coordinated  peak  and  center  bonding  sites. 


Monte-Carlo  calculation  of  the  distribution  of  adatoms 
on  a 30  X 30  lattice  for  Model  2 at  two  different  tempera- 
tures. a)  275°K,  order  parameter  = .972,  b)  325°K,  order 
parameter  = .804.  Each  atom  has  had  on  the  average  700 
chances  to  take  a step. 


MODEL  I 


€a=  Cgf  - 0.069  eV 
€[,=  + 0.15  eV 


• K MODEL  2 

K ea=-0.072eV 
K • eb=  + 0.080eV 

• ec=-0.049eV 


MODEL  3 
€a  • “ 0.09  e V 
€b*  €c*  + 0.075eV 
6^,*  0.03eV 
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